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wavelength. We show that a probe cylinder can still be cloaked at finite wavelengths, 
but the cloaking cylinder develops a non-zero scattering cross-section. We show that 
this latter effect is dominated by a monopole term in the case of an ideal (lossless) 
cloaking material, and by a dipole term in the case of a realistic (lossy) material. It 
can be reduced but not eliminated by variations of geometric or dielectric parameters 
of the cloaking cylinder. 
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1. Introduction 

There is much current interest in the possibility of cloaking or hiding objects from 
scrutiny by electromagnetic waves. At least three techniques have been proposed to 
achieve this: one avoids detection by surrounding the target body with a metamaterial 
shell which guides light around the central cavity [H [2], the second again relies on a 
metamaterial [3j HI El El [7J [8] , which this time cloaks by resonance an external region, 
while the third uses a structured metamaterial which provides cloaking by in effect 
folding space back upon itself [91 [10] . The technical challenges of making such systems 
in practice are enormous, but it should be realized that these and other proposals for 
cloaking offer complementary characteristics, which implies that work on a range of 
them is valuable. 



Table 1. (Jomp, 


arison or three cloakin 


g methods. 




Mechanism 


Refraction 


Reaction 


Unfolding 


Region 


Internal 


External 


External 


Structure 


Metamaterial 
Shell: e, n > 
vary with position 


Metamaterial 
homogeneous 

£s + £m = 


Metamaterial 
Shell: e, fj, > 
vary with position 


Equations 


2D, 3D Maxwell 


2D quasistatics 


2D, 3D quasistatics 


Experiment 


Yes 


No 


No 


Problems 


Bandwidth, 
structuring shell, 
energy dissipation. 


Bandwidth, 
achieving s s , 
scale size <C A. 


Bandwidth, 
structuring shell, 
scale size <C A. 



We present a number of salient characteristics of the three methods in Table 1. 
Cloaking by refraction requires a structured metamaterial shell to divert light around 
a cavity in which the object to be hidden is placed (internal cloaking). It is designed 
using a full solution of Maxwell equations in two or three dimensions, and there has been 
an experimental demonstration of this mechanism in the former case The second 
and third methods offer the complementary feature of concealing a body in a region 
close to, but outside, the cloaking system (external cloaking). In the case of cloaking 
by reaction, the object is concealed by virtue of a plasmonic resonance, which requires 
the material in the cylindrical cloaking shell to have a dielectric constant close to the 
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negative of the dielectric constants in the core and matrix regions surrounding it. It 
has been studied to this point mainly in two dimensions. The third method is the most 
recent, and has features in common with each of the other two. It achieves cloaking by 
using a spatially varying dielectric constant and magnetic permeability, designed using 
the same principles of transformation optics at a basis of refractive cloaking. The goal 
here, however, is to in effect fold space back upon itself, and, loosely speaking, to hide 
the object within the enfolding. 

It is our purpose here to study one problematic aspect of cloaking by reaction, which 
is implicit in two previous papers [7J [12] , but which has not been studied systematically. 
The problem does not arise for cloaking by refraction or unfolding, by virtue of their 
different mechanism of operation, which ensures that both the cloaking system and the 
cloaked object are hidden to an equal degree from electromagnetic probes. As we shall 
see, this is not necessarily for cloaking by reaction, where it is quite possible for a larger 
cloaking system to successfully cloak a small object, but to be itself quite visible. Using 
a visual analogy from the animal kingdom, we refer to this as " Ostrich Effect" : the large 
object hides the small object, but the large object does not hide itself. The possibility of 
the Ostrich Effect, was signaled in a paper by Milton, Brian and Willis p3]: "Besides 
invisibility there is what we call cloaking where the surrounding material does not have 
to be carefully adapted to suit the object to be made invisible. The cloaking device may 
be invisible or visible, although obviously the former is more interesting." 

We take the viewpoint here that for many purposes, the Ostrich Effect will be 
undesirable, and so we provide examples of the effect, explore its underlying physics, 
and quantify the circumstances under which it is greatly reduced. In Sec. II we present 
two figures taken from simulations showing the quenching of the dipole moment of a 
probe cylinder in the vicinity of a cylindrical shell with realistic values of the complex 
dielectric constant. These figures illustrate clearly the Ostrich Effect. In the next 
Section we consider a coated cylinder interacting with an incident plane wave and solve 
this scattering problem in closed form. We also take the long wavelengths limit of the 
formulation, in order to exhibit the transition from dynamics to qualitative. In Sec. IV 
we use the scattering cross section of the coated cylinder to show the counterintuitive 
result that a small imaginary part of e s actually benefits cloaking at finite wavelengths, 
since it makes dipole rather than monopole terms dominant in the scattering cross 
section. This is in keeping with the results of Hao-Yuan She et al [H] but not with 
those reported by Min Yan et al [T5] . 

2. Description of cloaking numerical simulations 

Let us consider a two-dimensional physical system comprising a coated cylinder centred 
about the origin of coordinates and a probe (solid) cylinder on the y-axis. Both cylinders 
are perpendicular to the xj/-plane. The shell and core radii of the coated cylinder are, 
respectively, r c = 20nm, r s = 65nm, while the radius of the probe cylinder is a = 5nm. 
Also, the core and shell relative permittivities are e c — 1, £ s = — 1 + 0.1 i. The relative 
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permittivity of the probe cylinder is e = e s and the relative permittivity of the matrix 
is £ m = 1. All the components are non- magnetic so that the relative permeabilities are 
/z c = /z s = // = // m = 1, where /x is the permeability of the probe cylinder. 

MO I I^^^m iiiiMi'iii'iiiiiii'ia^MMi'iMM ■imim im^m'i ' ' ' 




Figure 1. Left: Contour plot of \H Z \ as a function of position for a system consisting 
of a coated cylinder (r c = 20nm, r s = 65nm, e c — 1, e s = —1 + 0.1 i, /i c = fi s = 1) 
interacting with a probe cylinder (a = 5nm, e = —1 + 0.1 i, fx = 1), and irradiated 
by a H z polarized plane wave with wavelength 600nm coming from above. The probe 
cylinder is within the cloaking region bounded by the dashed circle, at a distance of 
85nm from the origin. Right: magnitude of the dipole moment of the probe cylinder 
as a function of its position indicated by the red line. The magnetic field varies in the 
range 0.63 < \H Z \ < 1.79, whereas the incident plane wave is normalized to \H z mc ^ | = 1 

This physical system is subjected to an incident plane wave having A = 600nm 
and with the wave vector in the xy-pl&ne (in-plane incidence) and polarized with the 
magnetic field parallel to the cylinder axes (H z polarization). 

The probe cylinder is polarisable, and has a dipole moment proportional to the 
total electric field at its position. When the probe cylinder moves along the y-axis and 
enters the cloaking region, marked by the dashed circle in figures [I] and [2j one can see 
the effect of cloaking, in the sense that the dipole moment of the probe cylinder tends to 
zero within the cloaking circle of radius r# = v^( r s/ r c) IS Ej (see figure [U right panel). 
Consequently, the probe cylinder is successfully cloaked within r # but not of course 
outside it, compare the left panels of figures []] and El Nevertheless, the cloaking system 
is not invisible, since the coated cylinder distorts the incident plane wave. Note that 
figures [H and [2] are frames in the animation available with this paper, which illustrates 
the quenching of the dipole moment of the probe cylinder within the cloaking circle, 
and its re-emergence outside it. 

The relative permittivity value chosen for the shell in figures [1] and [2] is comparable 
to that of silicon carbide near A = 10/im [16], and has an imaginary part somewhat 
lower than that of silver at the wavelength in the ultraviolet where the real part of its 
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Figure 2. As for Fig. Q] with the probe cylinder outside the cloaking region. Now, the 
probe is at a distance of 200nm from the origin, and 0.21 < \H Z \ < 1.81. 



permittivity passes through -1. While we refer to distances in nanometers, in practice the 
relevant parameter is the wavelength of the incident radiation divided by a characteristic 
length, say the outer radius r s of the coated cylinder. Hence, the results shown in figures 
[1] and [2] and subsequent figures can be easily applied to systems rescaled to correspond 
to materials other than those mentioned. 

Note that in the previous figures and simulations we have presented [6j [7J [8] 
to illustrate resonant cloaking, the value of the imaginary part of e s was chosen to 
correspond to the mathematical analysis, rather than to practical materials. We have 
found that good quenching of the dipole moment of the probe particle can be achieved 
even with quite significant imaginary parts for e s , provided r# is sufficiently in excess 
of r s , so that the probe particle can move deep within the cloaking region. 

Despite the quite effective cloaking of the probe cylinder shown in the right panels 
of figures [H and [2], the left panels illustrate strong distortion of the incident wave in 
the vicinity of the cloaking cylinder. Such variations of magnetic field strength would 
compromise any attempts to hide the compound system of cloaking cylinder plus probe. 

We mention that in all numerical computations we have used the scattering-matrix 
method pT], which is based on the expansion of the fields in terms of Fourier-Bessel 
series around each cylinder. By using the scattering matrices of each cylinder and the 
translation properties of Fourier-Bessel functions, the method leads to the inversion of 
a linear set of equations. 



3. The quasistatic limit 



For a coated cylinder centered at the origin of coordinates, we represent the electric 
and magnetic fields E z and H z (denoted here by V), by series expansions in terms of 
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cylindrical harmonics 



f A\J,{k c r)e l 



e=-oc 



'ATJi(k m r) + BrH^\k m r 



> e 



i(l3z~u}t) 



Jtt> 



where Jg(.) and H^\.) represent the Bessel and Hankel functions of the first kind. The 
three forms of the series expansions in (Tj[|) correspond to the domains < r < r c (inside 
the core of the coated cylinder), r c < r < r s (inside the shell of the coated cylinder) and 
r > r s (in the matrix), respectively. Also, the superscripts c, s and m label the fields 
inside the cylinder core, cylinder shell, and in the matrix, respectively. Thus, we have 
the wavenumbers k 2 = u 2 e c fi c , k 2 = u 2 e s fi s , and k 2 ^ = u 2 e m fi m . 

The function V has to satisfy the boundary conditions, i.e., the continuity of the 
tangential components of the electric (E z and Eg) and magnetic (H z and Hg) fields 
across the core and shell surfaces. When the coated cylinder is subjected to an incident 
radiation which is perpendicular to the axis of the cylinder, we have (3 = 0, and the 
problem can be reduced to solving two independent problems [T9] : 



• E z polarization, when H z = and the transverse parts of H are generated by V E z , 
and 

• H z polarization, when E z = and V H z gives the transverse components of E. 

In the present analysis we are interested in the relation between the coefficients in 
the matrix, which has the form 



Am 



-MpBT 



L e r> £ . (2) 

The coefficients A™ are determined by the sources of the field applied to the structure, 
and satisfy the field identity [8] 



oo 

E 



Af J e {k m r) e 1 " = source field . 



(3) 



Hence, we obtain the coefficients Af by expanding the source field in terms of cylindrical 
harmonics Je(k m r) e lie . 

Here, we also consider that the field applied to the physical structure is a plane 
wave field. In cylindrical coordinates, for H z polarization (H x = H y = 0), a magnetic 
plane wave is described by the formula 

where ipo is the angle of incidence with respect to the x-axis. We consider the case of 
in-plane incidence ik z = (3 = 0) so that the exponential in (jlj) can be expanded in terms 
of Bessel functions of the first kind 



(5) 
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Consequently, for a coated cylinder subjected to a plane wave incoming field, 
perpendicular to the cylinder axis, we have the coefficients 

A™ = H i n e- [n ^° , (6) 

where 

{7r, if the radiation comes from x = +00, 
(7) 
0, if the radiation comes from x = —00. 

In the case of E z polarization we obtain an equation identical to ([6]) for the 
coefficients of the electric field. 



3.1. H 7 Polarization 



We concentrate now on a coated cylinder, centered about the origin of coordinates, 
made from non-magnetic materials for which /i m = // s = // c = A'oj so that e m = n^e , 
e s = n^so and e c = nle , where e is the dielectric constant of free space, and n ; (i 
= m, s, c) represent the refractive indexes of the matrix, shell and core, respectively. 
The boundary conditions coefficients M% from (T5]), are derived by eliminating A\ in the 
equations 



/tm 



x 



Je(k m r s 
Z s J^k s r c ) 



H £ (k m r s 

He(k s r c ) 
Z s H' e (k s r c ) 



Je{k B r B 
Z s J f £(k s r s 



H e (k s r s 
Z s H' e (k s r s 



Je(k c r c ) H £ (k c r c 
Z c Jft(k c r c ) Z c Hg{k c r c 













H, 



are Bessel and Hankel functions of the first kind, the 



where J^(.) and Hi 

prime indicates the derivative of the corresponding function, and Z\ = \fji\je\ (i = 
m, s, c) represent the impedances of the matrix, shell and core, respectively. For E z 
polarization we obtain the relation between AJ 1 and BJ 1 by changing Z x — > 1/Zi in (JHJ). 

In the quasistatic limit (k m — > 0), we approximate the Bessel functions by the first 
term in their series expansion, i.e. 

■ Z \ 2 

.2, 



Hz) « 1 - (; 



Uz) 



1 /z\ n 
~n~\ V2, 

f-ir 



H (z) « J (z)+i- 



(-n)! V2 
2 



for n > 0, 



for n < 0, 



7T 



7 + log 



-- ( - ) (n - 1)! for n > 0, 



7T \ Z 



I (2 . (_ n _ -q, for n < Q) 
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7T yz y 

where 7 s is the Euler-Mascheroni constant [20]. Then, we substitute these expressions 
in (jSj) and write as a fraction. Using the limit k m — > we determine the coefficient of 
fc^ in the numerator and the coefficient of in the denominator. Thus, the quasistatic 
limit of M t for £ > is 

i / 2 \ 2 ' 

m, = m_^ — — eiie-iy.jt, (9) 

7T \ fc m r s J 

where 

_ (e s - e c )(e m - gg) + rf (e 8 + e c )(e m + e s ) 
£ r^(e s - £c ) (e m + e s ) + rf(e s + e c ) (e m - e fl ) ' 
For £ = 0, we obtain a completely different form 

M ° 7T (/c m r s ) (e s - e c )(r c /r s ) 4 + (e m - e s ) ' ^ 

Note that in all these calculations we made no assumption about the nature (real or 
complex) of permittivities or refractive indices. 

To relate the long wavelength limit of the dynamic problem with the corresponding 
problem in electrostatics we apply the same method as in Ref. [2T]. Thus, the boundary 
conditions for our problem correspond to an electrostatic problem in which the inverse 
of the dielectric constants (e c — > l/e c , e s — > l/e s and e m — > l/e m ) have to be considered. 
This will also change 7^ — > —7^. Now, the boundary conditions for £ 7^ can be 
written in the form 

2C 



— ) (12) 



Note that, here, we separated the product k m r s , which is dimensionless, so that k m and 
r s are considered as multiplied, respectively divided, by a length unit. 

In electrostatics, the corresponding relationship between the coefficients Ai and Be 
which controls the response of a coated cylinder to an external field, has the form 0, H] 

At = li%. (13) 

'"s 

Now, by comparing (|12p with f[T5j) we may infer the relation between static and dynamic 
multipole coefficients 

S < * if) 21 W - ^ 5 ™ * TTTT ^ > for ^ °' ( 14 ) 
n \k m J Ji{k m ) 

Note that fc m is dimensionless according to the note after (1121) . 
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In electrostatics, the partial resonances of a three-phase composite consisting of 
coated cylinders are defined by the equations (3J H] 

e c + e s =0 (core — shell resonance), (15) 

£ s + £ m = (shell — matrix resonance), (16) 

when ffl3l) becomes 

A, = e ±±±*, ( 17 ) 

or 

Ai -e m -s c (18) 
respectively. In the first case, the field inside the coated cylinder is exactly the same 
as would be found within a solid cylinder of radius r s and dielectric constant e c , while 
the potential outside the coated cylinder, in the matrix, is precisely the same as that 
outside the solid cylinder [SI El- The second case corresponds to a solid cylinder of radius 
r s/ r c > r s (the geometrical image of the core boundary with respect to the shell outer 
boundary), and dielectric constant e c . Now, the field external to the coated cylinder 
and beyond the radius r* = r^/r c is the same as that external to the solid cylinder [SIE]. 

Since it is the relationship between the coefficients Af 1 and Bf 1 which controls the 
response of a coated cylinder to an external field, equations (fl3l) and (fT2|) show that this 
response is determined by je in electrostatics as well as in the long wavelengths limit 
of electrodynamics. The limiting process is smooth and therefore, we expect a resonant 
behavior accompanied by cloaking effects, even for nonzero frequencies, when one of the 
conditions ffT5l) or ffT6l) is satisfied. 



3.2. E z Polarization 

Now, in the long wavelength limit, the boundary conditions coefficients Mi from ([2]) 
take the form 

M , = M _, „ _i ( JL ) 2M ™V±£5 > , (19, 

7T V k mr s J {e s - e c ) {r c /r s ) M+2 + (e m - e s ) 

for i ^ 0, and 

i 



M ~ — I - — ) ^ r. (20) 



Note that there are no terms of the form e c + e s or e s + e m in (1T9"]) or ( f20l) to indicate 
a core-shell or shell-matrix partial resonance. Again, the limiting process is smooth 
and, consequently, we do not expect a resonant behaviour of the coated cylinder, for 
any frequency, in the case of E z polarization. Consequently, when the coated cylinder 
is irradiated with a field of a general polarization, that is a mixture of H z and E z 
polarizations, or in the case of conical incidence, the cloaking by resonance will never be 
perfect, or even may be completely ruined, due to the contribution of the E z polarized 
component. 
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Resonances similar to those in Sec. [S] can occur in the case of a coated cylinder 
made from magnetic metamaterials with permitivitty Eq and permeabilities fi c , /i s , and 
/i m . Now, the boundary conditions coefficients Me from (j2J) take the forms ([9]) and 
( TTTj) with £i replaced by fa [18j [H]. Hence, the magnetic partial resonances of coated 
cylinders are defined by the equations 

fa + A*s =0 (core — shell resonance), (21) 
A* s + A*m = (shell — matrix resonance). (22) 



4. Attempts to minimize the Ostrich Effect 

As a measure of effectiveness of cloaking we choose the total scattering cross section. 
For two-dimensional problems, the total scattering cross section is defined as the ratio 
of the total power scattered by an object, to the incident power per unit length [T9l 122] 

. oo . oo 

^ r E i^ m i 2 = r |5 ° m|2 + '^ n|2 • (23) 

"Til a "'m ™m e -. 

£=— oo 4=1 

Here, the Bf 1 coefficients have the exact form given by (jHJ), that is 

B? = ^ Af, (24) 

where 

P = {[Z m Jt(k m r s ) Je(k s r s ) 
x [Z c Hi(k s r c ) J' e (k c r c ) 
- {[Z s H' e (k s r s ) Je(k m r s ) ■ 
x [Z s J' e (k a r c )J e (k c r c ) 



- Z s J e (k m r s ) J e (k s r s )} 
-Z s H' e (k s r c )J e (k c r c )}} 

H e (k s r s ) J' e (k m r s )] 

- Z c JAk s r c ) J'Jk c r c )]} , 



Q = {[Z m H' e (k m r s ) J e (k s r s ) - Z s H t {k m r s ) J' e (k s r s )] 
x [Z s H' e (k s r c ) Jt{k c r c ) - Z c H e (k s r c ) J' e {k c r c ))} 
- {[Z m H' e (k m r s ) H e (k s r s ) - Z s H t {k m r s ) H' e (k s r s )} 
x [Z B J' e (k s r c ) Je{k c r c ) - Z c J £ (k s r c ) J' e (k c r c )]} , 

and with A™ from ([6]) for E = 1 and if) = 0. 

In the case of the resonance e c = e m = 1 and e s = —1, numerical simulations show 
that by using the form (1241 and the series fl23|) truncated to iV trunc = 6, we have 

af— ) = ± "g° |5n 2 - f |5 m ! 2 = ^, (25) 

*- — iV trunc 

starting at about A = 10 r s (see figure[3]). From the same wavelength up, the contribution 
of the dipole terms, given by cr^ — o~[°\ becomes very small. 




If we use the expression of B™ in the quasistatic limit (see Sec. 13. ip . with A™ from 
E = 1, and ip = 0, we obtain 



B) 



1 7T 



e c ){r c /r s Y + (e x 



A 



1 



2/' 



for 



for 



0. 



> 1. 



(^-1)!7£ 

In the case of core-shell-matrix resonance, that is e c + e s = and e s + e m 
coefficient 7^ defined in f fTUj) tends to infinity so that, for i > 1, — » 0. 
situation arises when 



(26) 

= 0, the 
Such a 

(27) 



fas in the case of e r 



1 and e s = —1), and the total scattering cross section is 



determined only by the zeroth-order multipole 



a 



QS 



4 IR' 

F 1 



(2f 



which tends rapidly to zero as the wavelength increases. Actually, B™ from ( 12^1) tends 
very slowly to the form (1261) . This last form has been obtained by taking the first term in 
the series of all Bessel functions, except Jq(z) and H (z)). For complicated expressions 
like ([21]) the series expansions require more terms for accuracy, as they contain products 
of four Bessel functions. 

The main result is that for long wavelengths a t is determined by the 5™, only. 
The dominance of zeroth-order multipole is also present in the case of coordinate 
transformation method [TJ. This case has been analyzed by Yan et al [T5] . 

Here, we have considered that e s = —1 is real, which is unphysical. Physical 
materials with negative permittivity (usually metals) are lossy so that we have to 
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consider a coated cylinder with the shell material having a complex permittivity 
e s = — 1 + id, where 5 > determines the loss in the shell. A detailed analysis of 
this case shows that, in the limit of long wavelengths, the total cross section at is now 
dominated by the dipole coefficients B± v 




Figure 4. Total cross section a\ (blue curve), a\ (red curve) from ([^5)) , and 
(Tj 1 ^ — (green curve) as functions of the wavelength of the incident plane wave, for 
a coated cylinder (r c = 20nm, r s — 65nm) at resonance (e c = e m = 1, £ s = —1 + 0.1 i, 
Mc = Ms = Mm = 1 )■ 



We start with the analytic form of B™ coefficients ( I24J) and set the factor A™ = 1. 
In fact, A™ defined in has the modulus = Eq and we consider Eq = 1. Now, for 
a finite S, from the exact form (|24|) we obtain the following series expansion 

5 m = a 4 {6) k* + a e (S) k & + O(k 8 ) , (29) 

where 

a 4 (5) = i^(r c 4 - r s 4 ) + ^(r c 4 - r s 4 )5 , (30) 
a 6 (5) = i^-(r c 6 - 2r c 4 r s 2 + r s 6 ) + -^-(7r c 6 - 12r c 4 r s 2 + 5r s 6 )5 - i-^-(2r c 6 - 3r 4 r s 2 + rf )<5 2 . 



(31) 



32 v c a , lg2 v c C B ■ B , 

It is easy to check that the first term in the limit 

fim B™ = i ^ (r c 4 - r s 4 ) A; 4 + i ^ (r c 6 - 2r c 4 r 2 + r s 6 ) A; 6 + O(k 8 ) , (32) 

is of the form (1261) for £ = 0, if we set e c = £ m = 1 and £ s = —1 in (T26|) . Consequently, 
the behaviour of I?™ as function of k Q and 5, in the domain of long wavelengths, can be 
summarized as 

f a 4 (5)fc 4 if 5^0, 

Bo * (33) 
[ a 4 (0)fc 4 if 5 = 0. 

Now, we analyze the dipole term B™. Firstly, from (1241) we obtain 

B? = b 2 (8) kl + h(S) fc 4 + O{k G ) , (34) 
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where 

h(s ,_* rl{rl-rl){2-i5) 
2[ ) 4 -r 2 5 2 + r 2 (2 + i5) 



7rrs 2 (r 2 -r 2 ) vr r 2 

7T~^ 0+1 



■ r 



16 rl 



[r 2 - r 2 )6 2 + 0(S 3 



vr 4 . / r, 



6 4 (5) = -i-r s 4 log ( ^ 



?r s 4 log(^)rt 



We also have 



limR m 

5^0 



(35) 
(36) 

(37) 



Finally, the behaviour of Bf 1 as function of ko and 5, in the domain of long wavelengths, 
can be summarized as 



b 2 {5)k 2 if 5^0, 
64(0)^ if 5 = 0. 



(38) 



It follows that, in the long wavelength limit, the scattering cross section is 
dominated by the monopole term = 4 \B™\ 2 /k m if 5 = 0, and by the dipole term 



— <j [ t v> = 8 \B^\ 2 /k m when 5^0. This last result agrees with that obtained by Alu 
and Engheta [5] who have also shown that the dipole term dominates the scattering 
cross section, for lossy materials. 





Figure 5. Scattering cross section a[ 6 \ from equation (|23|) . as a function of wavelength 
(A) and core refractive index (n c ). Left: n c < 1. Right: n c > 1. The parameters of the 
physical system are: r c = 20nm, r s = 65nm, e c = 1, e s = — 1 + 0.1 i, /i c = /i s = 1- Red 
indicates large cross sections (up to 18.0), while blue indicates smaller values (down to 
0.13). 



Figure H] shows the cross section a t as a function of wavelength calculated now for a 
realistic value e s = — 1 + O.li of the dielectric constant of the shell. In comparison with 
figure El we see that the dominant contribution to the cross section now comes from the 




dipole terms rather than the monopole terms. Despite this difference, the cross section 
is well approximated by its leading term when the wavelength reaches around 10 times 
the shell radius. Note that in figure [3] the cross section varies as approximately 1/ A 2 in 
the region of A between 600 and 1000; this is far from the quasistatic behaviour of 1/A 7 
expected from equation ( 1331) . showing that for this ideal case the cross section is not 
well represented by quasistatics even at A = 1000. By contrast, for figure 0] the cross 
section goes as 1/A 3 , in line with the quasistatic estimate (see equation (1381) . and also 
Panicky and Phillips [T9]). 

We can examine whether the Ostrich Effect can be reduced by making appropriate 
choices of the free parameters of the cloaking system: n c , r c and r s . We study the effect 
on the scattering cross section of varying these parameters in figures [5] and El Figure 
|5] shows the effect of varying the core index, both below and above the value of unity 
used in previous figures. While the effect of n c varies with wavelength, in general one 
sees from figure [5] that values around unity deliver the lowest cross sections. In figure [6] 
we study the effect on cross section of varying radii. Here, the cross section values have 
been normalized by dividing by the cylinder diameter, to give a dimensionless value. 
The geometric parameters r c is shown divided by the wavelength, while the horizontal 
axis gives the cloaking radius r# divided by r s . The leftmost curve gives the contour 
on which the cross section is equal to 1% of the geometric value. Along this contour, 
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if we want to have say r# = 2.5r s , to give a relatively large cloaked region, we need 
r c ~ 0.003A and r s ~ 6.25r c , or r s ~ 0.019A. These relatively strict tolerances illustrate 
the difficulty of achieving low cross section values at finite wavelengths. Figures [7] (a) 
and (b), show two field distributions corresponding to the probe inside and outside the 
cloaking region, with Im(e s ) now set to 0.01, in order to make more effective the cloaking 
action. We can see in these figures that we have achieved a satisfactory combination 
of effective cloaking of the probe and virtual elimination of the Ostrich Effect. We 
can quantify this by introducing a dimensionless quantity we call visibility defined in a 
similar fashion to the quantity in interference optics: 

v = (|# z {max) | - |# z (min) |) / (|# z (max) | + |# z (min) |) . 



(39) 



Here, |ifi mm ^| and |iY z vmdJt; | denote the minima and the maxima of modulus of H z in 
the region outside a circle of radius r#, if the probe is within the cloaking region, and 
outside the minimal circle containing the coated cylinder and the probe, if the probe is 
outside the cloaking region. The values of v are 0.526 for figure [3(a) and 0.021 for figure 
E^b), a satisfactorily small value. 



r (max) I 



< 

— 1 — Ss 


■ i i i i i i i i i i • i ■ 

> 




v. ; 



a 



(b) 



Figure 7. Contour plot of \H Z \ as a function of position for a system consisting of a 
coated cylinder (r c = 1.8nm, r s = 11.25nm, e c = 1, s 3 = —1 + 0.01 i, fi c — Ms = 1) 
interacting with a probe cylinder (a = 2nm, e = — 1 + 0.01 i, \x = 1), and irradiated by 
a H z polarized plane wave coming from above, (a): The probe cylinder is outside the 
cloaking region bounded by the dashed circle, at a distance of 150nm from the origin, 
(b): The probe cylinder is outside the cloaking region at a distance of 13.5nm from 
the origin. The visibility ([59)) takes the values v = 0.526 (a) and v = 0.021 (b). 
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5. Conclusions 

We have presented numerical results displaying clearly the tendency for it to be more 
difficult to hide the larger cloaking system than the smaller object it is trying to conceal 
from electromagnetic probing, and we have analyzed this effect to quantify to what 
extent it can be overcome. Our results are conveniently summarized in figure [7J, and 
give size limits on the cloaking system in terms of the wavelength. These size limits in 
fact just require both the cloaking system and the system it is cloaking to be in the 
quasistatic regime. 

We have also shown that this regime in fact sets in at shorter wavelengths for a 
resonant cloaking system with a small amount of loss, compared with the case of no 
loss. We have confined our studies to cloaking systems which have spatially uniform 
shells, but it may be the case that structured systems of the sort described by Farhat et 
al [23] may be designed which inhibit multipole responses in such a way as to ensure the 
onset of quasistatic behaviour at shorter wavelengths than indicated by figure [7J This 
would be valuable in possibly simplifying the construction of cloaking systems which 
operate by plasmonic resonance, while it would be bringing their geometry closer to 
that of systems which cloak by refraction. 
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